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Abstract 

Evidence for free precession has been observed in the radio signature of several pulsars. Freely 
(N ■ 

precessing pulsars radiate gravitationally at frequencies near the rotation rate and twice the rotation 



rate, which for rotation frequencies greater than ~ 10 Hz is in the LIGO band. In older work, the 
gravitational wave spectrum of a precessing neutron star has been evaluated to first order in a small 
\ precession angle. Here we calculate the contributions to second order in the wobble angle, and we 



find that a new spectral line emerges. We show that for reasonable wobble angles, the second-order 



line may well be observable with the proposed advanced LIGO detector for precessing neutron stars 
as far away as the galactic center. Observation of the full second-order spectrum permits a direct 

•i-H . 

measurement of the star's wobble angle, oblateness, and deviation from axisymmetry, with the 



potential to significantly increase our understanding of neutron star structure. 
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I. INTRODUCTION 



For some decades now, freely precessing neutron stars have been considered as possible 
sources of detectable gravitational radiation. Early treatments include the ones by Zim- 
mermann et al. and by Alpar and Pines [4J. More recently interest in the subject 

was revived by compelling observational evidence for free precession in the radio signature 
of PSR B1828-11 fl as well as PSR B1642-03 Q. Evidence for free precession had also 



been put forward for PSR B0531+21 (the Crab pulsar) (ft |8| and PSR B0833-45 (the Vela 
pulsar) jj. These examples suggest the existence of more such objects, some of which might 
be viable sources for current or future ground-based gravitational wave detectors such as 
LIGO. 

Prompted by the evidence for free precession in PSR B1828-11, Cutler and Jones 
revisited earlier work which seemed to suggest that when a neutron star is treated as an elas- 
tic body, gravitational wave back-reaction will quickly damp the precessing motion. Careful 
consideration led them to conclude instead that in realistic precessing neutron stars, grav- 
itational wave damping will play but a minor role in the star's evolution, although there 
will be other sources of da^in, Jones aad Anders Q then _ t ed a detailed 
model of freely precessing neutron stars and compared it with the observational data. They 
concluded that the precessional interpretation of these observations required superfluid vor- 
tex pinning of the star's crust to its core to be inhibited, thereby challenging some of the 
then-prevalent ideas concerning a neutron star's liquid interior. Nevertheless, an analysis by 
Link and Epstein supported the precessional hypothesis for PSR B1828-11. 

In jl3|, Jones and Andersson evaluated the viability of precessing neutron stars as sources 
for gravitational wave detection in the near future. To estimate the strength of the signal, 
they used the waveforms computed by Zimmermann et al. . Their conclusions appear 
to imply that gravitational waves from precessing neutron stars are unlikely to be detectable 
even with an advanced LIGO detector. However, their paper was written before the more 
recent analysis of the expected capabilities of LIGO II by Fritschel [14I which is the one 
we will use. 

Here we re-examine the gravitational waves generated by rapidly rotating, precessing neu- 
tron stars, evaluating whether the radiation from such sources might indeed be detectable, 
and what we may learn from their detection. 
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The motion of a freely precessing neutron star has two characteristic frequencies asso- 
ciated with it: the precession frequency and the rotation frequency. The corresponding 
gravitational wave spectrum consists of a series of discrete spectral lines, with a princi- 
pal line at or near a harmonic of the rotation frequency, surrounded by companion lines 
separated tan the main line by .nn.tip.es of the preeessiona. tajnenc, 

In |3J , Zimmermann evaluated to first order in small wobble angle the quadrupole gravita- 
tional radiation emitted by freely precessing, non-axisymmetric neutron stars, approximated 
as rigid bodies. At first order the spectrum includes two spectral lines: one at approximately 
the pulsar frequency and another at its first harmonic. Unfortunately, knowledge of the 
first-order spectrum yields little or no information about the physical characteristics of the 
neutron star emitting the radiation. From an observational point of view, this makes it 
imperative to study the spectrum in more detail. 

Our goal here is threefold. First, we want to show that it is reasonable to take the 
precession angle as the primary expansion parameter among the various 'small' parameters 
that characterize the problem. Secondly, we want to determine the conditions under which 
the second-order contributions to the gravitational wave spectrum of a precessing neutron 
star are observable. Finally, we would like to explore what observation of these second-order 
corrections can reveal about the physical properties of the source. 

Throughout this paper, we will use units where G — c — 1 unless stated otherwise. 
Wherever estimates are needed for the mass and radius of a neutron star, we will use 
M = 1AM & and R = 10 6 cm. The neutron stars of interest to us are fast-spinning ones; 
in estimates involving rotation frequencies we use 500 Hz as a reference. We shall treat 
neutron stars as rigid bodies and limit ourselves to the quadrupole approximation to the 



3 



16 



171 ]. i.e., we 



7j. In section |TT] 



gravitational radiation. When describing rigid body motion we will follow 
use the physical conventions of and the mathematical notation of 0, 
we review the closed expressions for the gravitational waveforms found in [| and in IIHI the 
first-order approximation of [3j. In section HV| the closed expressions are used to construct 
a second-order expansion, and a new contribution to the gravitational radiation spectrum 
of precessing neutron stars is found. Section IS first discusses the detectability of all three 
spectral lines. We then explain how physical information can be extracted from the second- 
order spectrum. Finally, section IVII provides a summary of the results and conclusions. 
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II. GRAVITATIONAL RADIATION AND PRECESSION 



The motion of an isolated, freely precessing rigid body is most conveniently described 
in terms of the Euler angles 9, <fr and ip that orient the body with respect to an inertial 
coordinate system. Following the conventions of denote the coordinates of the inertial 
coordinate system centered on the body by X, Y, and Z, with the coordinate vector ez 
in the direction of the body's angular momentum. Similarly, denote the coordinates of 
the body-fixed frame as x, y, and z, with coordinate vectors e x , e y and e z parallel to the 
eigenvectors of the body's moment of inertia tensor, and with the corresponding principal 
moments of inertia satisfying I z > I y > I x . The Euler angle 9 is then the angle between 
e z and e z , <fi is the angle between e x and the line of nodes that is the intersection of the 
ex x ey plane and the e x x e y plane, while ip is the angle in the e x x e y plane between the 
line of nodes and e x (cf. [Iy, figure 47]). 

Denote the time-dependent rotation matrix that transforms from the body-frame to the 
inertial frame by Rj^, where Latin indices refer to the inertial frame and Greek indices refer 
to the body-fixed frame. The quadrupole gravitational radiation from the precessing body 
is proportional to the second time derivative of the moment of inertia tensor in the inertial 
frame, which depends on the principal moments of inertia and the time-dependent rotation 
matrix. 

The standard quadrupole formula for the transverse-traceless (TT) gauge metric pertur- 
bation is 

where 

i jk = J d 3 X (xjX k - l -^ 2 5 3 ^j p(X), (2.2) 

with p(X) the density of the source. Eq. ()2.1|) can be written in terms of the rotation matrix 
Rjn, the components f2 M of the angular velocity in the body-fixed frame, and the eigenvalues 
of the rigid body's moment of inertia tensor Ji, J 2 , and /3 jjj: 

TT 2 T 

= ~R"j\*-^v-vR'vki (2-3) 

where summation over /i and v is understood. A^ v is given by 

A u = 2(A 2 Q 2 2 - A 3 Q 2 3 ), (2.4a) 
A 12 = TsQ^, (2.4b) 



with the other components being defined by symmetry and cyclic permutation of the indices. 
A^ and are shorthand for 



Ai = h ~ h, (2.5a) 
r x = A 2 - A 3 + ^i, (2.5b) 

where the other components are again given by cyclic permutation of the indices. 

To facilitate the description of the radiation, introduce two orthogonal unit vectors v and 
w with vxw directed toward a distant observer of the emitted gravitational waves. Without 
loss of generality we may assume that the observer lies in the ey x ez plane. Let i denote 
the angle between ez and the line of sight to the observer, and define 

v = ey cos(z) — ez sm(i), (2.6a) 
w = -e x . (2.6b) 

In this frame the TT gauge metric perturbation corresponding to the observed radiation can 
be described in terms of two scalar functions, h + (t) and h x (t), representing the radiation in 
the two gravitational wave polarizations: 

hj?(t) = h+(t) (e + )^ + h x (t) (e x ) y , (2.7a) 

with 

e + = v®v — w®w, (2.7b) 

e x = v <8> w + w g) v , (2-7c) 

K = \h^ + k , (2.7d) 

h x = l -h^e^. (2.7e) 

Using equation ()2.3|) for the TT gauge metric perturbation we have 

h+ = -- \(R y ^ cos(z) - R Zfl sin(z)) (R yu cos(z) - R zu sin(z)) - R XfM R xu ] A^ u , (2.8a) 
2 

h x = - (R yiJ ,cos(i) - R zlx sin(i)) R^A^. (2.8b) 

In these expressions time evolution arises from the time dependence of the angular velocity 
in the body frame and the rotation matrix. 
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When J 2 > 2EI2, with J the magnitude of the angular momentum and E the energy 



in the rigid body motion 
given by Q]: 



22j . the components of the angular velocity in the body frame are 



^3 



a cn(r, m), 

\ h (h ~ h) 
a [h (h ~ h) 
6dn(r, m), 



1/2 



where r is a rescaled time variable, 

r = bt 



(I 3 -J 2 )(J 3 -/0 
hh 



(2.9a) 

(2.9b) 
(2.9c) 

(2.9d) 
(2.9e) 



1/2 



cn, sn, and dn are elliptic functions [17|], and m is given by 

(J 2 - J^a 2 
(J 3 - h)hb 2 

The constants a and 6 are such that at t — 0, Oj = a and fi 3 = 6. When the precession 
angle is small, it will approximately equal a/6. 

The rotation matrix Rj^ is given in terms of the Euler angles 9, <p, ip. The time dependence 
of these can in turn be expressed in closed form in terms of elliptic and theta functions: 

hb 



cos 9 



J -dnr, 
h{h~h] 

Hh - h] 

'1 + 02, 



1/2 



cn r, m 



sn(r, m) 



with 



exp [2i(f>i(t)} 



04 - ™) 

2vrf 



(^2 

2tt 
T 7 



h 



2i d' 4 (ma) 



T 4 (i7ra) 

The constants a and T that appear above are given by 

il 3 b 



sn [2iaK(m)] 
T 



AK(m) 



hh 



(h-h)(h-h)\ 



1/2 



(2.10a) 

(2.10b) 
(2.10c) 



(2.10d) 

(2.10e) 
(2.10f) 

(2.10g) 
(2.10h) 
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where $4 is the 4 th Jacobi theta function with nome 



9 = exp l-^%^], ("DO 



and K(m) is the complete elliptic integral of the first kind 17]. 

Within the weak-field, quadrupole approximation these closed expressions fully describe 
the gravitational wave signal from a freely precessing, rigid body. However, after all is said 
and done, the problem involves just two periodicities determined by T and T'. Accordingly, 
the radiation should be resolvable into a discrete spectrum. Knowledge of the amplitude 
and frequencies of the spectral lines greatly simplifies the analysis of gravitational wave 
data. The discrete nature of the spectrum becomes apparent when h + and h x , as given in 
equations (|2.8jh are expanded in terms of 'small' parameters: the 'wobble angle parameter' 
a/b, the rigid body's oblateness, and its deviation from axisymmetry. Finding the second- 
order contributions to the spectrum in these parameters, evaluating the observability of the 
associated spectral lines, and exploring what can be learned from such observations will be 
our main considerations. 



III. FIRST ORDER CONTRIBUTIONS TO THE GRAVITATIONAL WAVE 
SPECTRUM 

In his first-order expansion of h + and h x , Zimmermann identified two expansion param- 
eters, m and 5: 



m = T F~T~ (t) ' (3.1a) 



' = • (31b) 

To find the 'first-order' contributions to the radiation spectrum, Zimmermann 

• expanded the rotation matrix R ifJL to first order in m and 5; 

• substituted the expanded rotation matrix into the relation between the inertia tensor 
in the inertial frame and that in the body frame; 

• substituted the inertia tensor in the inertial frame into the standard quadrupole for- 
mulae for the strains (J2.ll) ; 



7 



truncated the resulting h + and h x to first order in small a/b, small oblateness (I3 
h)/h, and small non-axisymmetry (J 2 — h)/(h — h)- 



This leads to 

h z = jmL (h-h) ^ + cog2( ^ cos(2firott) 

r i3 

7 3 (q ot + ^ pre c) 2 2/3 - (Ji + I 2 ) (ha\ . 
H jtt YT sm ( 2i ) cos[(i2 rot + U vrec )t\, (3.2a) 

r i 3 

2/3 (Oct + "prec)' 2/3 - (/1 + h) (ha\ • /.x • r/ n ^ n wi ,o 9M 
H £7 YT an ( Z ) Sm i( fi «rt + "prec^J- (3.2b) 

5" -^a \hv J 

The fundamental angular frequencies in the problem are 

2tt 



"prec — ^ 



TTfo f{h-I 2 ){h-h) 



2K(m) V Ii/ a 
2tt 2tt 

rot ^j-j^ ^j-j 



1/2 



(3.3a) 



j _( l + i*g*ay )IW _ (8JSb) 



Ji \ 7r^4(i7ra) 

The expressions (13.2)) represent Fourier expansions of the waveforms, with gravitational 
spectral lines at the angular frequencies 

n 1 = 20ot (line I), 

fi 77 = fi rot + fi pr ec (line II). 

Line I is a consequence of the body's deviation from axisymmetry. Non-axisymmetry is 
often discussed in terms of an 'equatorial eccentricity': if the body is approximated as 
an ellipsoid rotating about an axis associated with a principal moment of inertia, then it 
will look identical after half a rotation period, hence the frequency 2f2 rot . Line II encodes 
the precessing motion; recall that a/b is approximately the wobble angle. The coupling of 
precession and rotation causes the frequency of this line to be shifted upwards by f2 prec ; note 
that both motions are in the same sense. As we shall see later on, this coupling also has an 
effect near line I, although it only becomes visible at higher orders in the precession angle. 
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IV. SECOND ORDER CONTRIBUTIONS TO THE SPECTRUM 



Zimmermann's calculation of the first order contribution to the gravitational radiation 
from rigid body free precession [sj involved the expansion parameters m and 5 defined by 
the somewhat awkward expressions (|3.1a|) and (j3.1b|) . and at the end further truncations 
were made in terms of different small quantities. Here we will expand in terms of small 
parameters that have a more direct physical meaning. 

The wobble angle is only one of several small parameters in the problem. To set up a 
second-order expansion, we will first specify a suitable set of expansion parameters and then 
estimate the size of these parameters for realistic neutron stars. These estimates will tell 
us to what order we need to expand in parameters other than the wobble angle parameter 
a/6 to arrive at a sufficiently good approximation for the waveforms. They will also yield 
information about the amplitudes of the different spectral lines at first and second order. 

Unlike Zimmermann, who used the general expressions for the quadrupole radiation (j2.1|) . 
our starting point will be the closed forms ()2.8)1 for h + and h x . There is no significant 
difference between the two approaches, but expansions of the closed formulae will lead to 
somewhat more elegant expressions for the spectral line amplitudes. Overall, our aim is to 
make the derivation and physical meaning of both the first-order and second-order results 
as transparent as possible. 



A. Expansion parameters 

The parameter m defined by eq. f!3.1a|) encodes two pieces of information at once: the 
deviation from axisymmetry and the precession angle. To understand the role played by 
each, factor m as 

m = 16k 7 2 , (4.1) 
where the new small parameters 7 and k are defined by 

lhh ~ h (4.2b) 



16 h h ~ h 

The parameter 7 describes the precession angle while k describes the non-axisymmetry 
(i~ 2 — h) relative to the axisymmetric non-sphericity (I3 — 12). Note that in all cases that 



are of astrophysical interest, J1//3 — 1; this and the factor of 1/16 in k have been included 
for computational convenience. The parameter 5 (eq. 13. lb)) is related to the new parameter 
k; it will be discussed in subsection C below. 

We would like to stress that specifying these different expansion parameters from the 
outset does not imply a departure from the approach of [3j|. Close inspection of that analysis 
shows that 7 and k are in fact used implicitly to truncate to first order in a/b at the end of 
the computation. Introducing them from the beginning merely facilitates the calculation. 

In order to develop appropriate parameters for the second order expansion, we need to 
find out what astrophysics is involved in the new parameters 7 and k, and any other relevant 
'small' quantities. Note that in realistic cases, 7 and k will not be independent, since the 
precession is determined for the most part by the loss of axisymmetry. Ideally, it would be 
a simple exercise in classical mechanics to find the relationship: one could start from a non- 
precessing, axisymmetric body and then study what kind of precessing motions are induced 
by various continuous deformations. However, we may assume that during a starquake other 
factors will be at play, such as torques induced by partial crustal pinning of vortices in the 
interior. To estimate their effect we would have to take recourse to the existing models of 
pinning, which seem to be contradicted by the observations of apparently free precession in 

reason, we will mostly disregard any a priori relation between 7 
and k and keep the discussion as general as possible. 

B. Astrophysics of precessing neutron stars 

Let us first look into some estimates for oblateness, non-axisymmetry, and the precession 
angle. 

Oblateness was described by Baym and Pines jlq] as 

^rM^o. (4.3) 

n 

Here (3 is the rigidity parameter, which depends on the equation of state; in [19| it was 
derived to be 

(3 ~ 1.6 x 1(T 5 . (4.4) 

For very young neutron stars, the zero strain oblateness €q will depend on the shape of 
the star when the crust first solidified; subsequent starquakes as well as less violent plastic 
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deformations will change it. We will take it to be of the order of the ratio of rotational 
kinetic energy to gravitational binding energy: 

Q 2 R 3 

where / ro t = ^vot/2n. We have set M = 1.4M and R = 10 6 cm, which are the values for 
mass and radius we will use throughout this paper. 
Eqns. (|Oj) . (JOJ), and (|4~5j) lead to an oblateness of 

h ~ h - 7.7X10- 7 f-%-V f^V (4.6) 



i 3 V500 Hzy Vio -5 , 

The implied rang e ag rees with estimates by Alpar and Pines |4j; see also the discussion in 
Cutler and Jones jlO(. 

To arrive at an estimate of k, we need information as to how a neutron star loses its 
axisymmetry; here we will follow the analysis by Link, Franco, and Epstein [20]. According 
to their model, a starquake is expected to induce fault lines at angles of 30 to 45 degrees 
relative to the equator. Crust material shifts along the faults and accumulates in mountains 
with a typical height of 10 -2 cm. We may assume 

I 2 - h ~ M m R 2 , (4.7) 

with M m the mass contained in a mountain and R the radius of the neutron star. One has 

Apr 

where AR is the change in equatorial radius due to the decrease in oblateness, L is the fault 
length, and M c is the mass of the crust. The factor in front of M c in ()4.8|) is the fraction of 
the star's area that moves with the fault. Assuming an outer crust density 23] of p ~ 10 10 

am 3 and a crust thickness in the order of R/10, we can use the expression for AR/R in 
to arrive at 

^^■^^K^^^te)" 1 ' (4 - 9) 

where t q is the average time between large starquakes and t age is the spin-down age. 
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We now have sufficient information to estimate k. Using (j4.2b|) . ()4.6|) . and ()4.9|) . we find 



k~2.2x10- 4 - ^— MrM • (4-10) 

The most important 'small' parameter in our analysis will be the wobble angle 7. For 
sufficiently small rotation frequencies, the wobble angle is essentially unconstrained. For 
rapidly rotating neutron stars (which are the ones of interest for gravitational wave detec- 
tion), Jones and Andersson ll( propose the following expression for the maximum wobble 
angle that can be supported by the crust: 

500 Hz\ 2 /M br eak\ 



7m „^, 8xl0 - 2 ^_j { ™ } , (4.11) 

where Wbreak is the strain at which the crust fractures. It is easy to understand this expression 
qualitatively; if f TOt is large, so will be the oblateness and hence the equatorial bulge. The 
star will then have to displace more matter while precessing, leading to higher crustal strains. 
In j^jj, the breaking strains of terrestrial materials were extrapolated to yield estimates in 
the range 10~ 4 < Mbroak < 10~ 2 , although the reliability of such extrapolations is unclear. 
Note that if the breaking strain is near the upper limit of this range, Wbreak — 10~ 2 , the 
expression (|4.11|) implies an unconstrained precession angle even at a rotation frequency of 
100 Hz. 



C. Setting up the second order expansion 

From the estimates ()4.(jj) . ()4.9j) . ()4.1()j) and ()4.11|) of the previous subsection, we infer that 
in the frequency band of interest, 



max 

(4.12) 



for /i, z/, p = 1, 2, 3. 

We will choose 7 as our primary expansion parameter. The approximation ()3.2|) for the 
waveforms is of first order in 7; here we will go to order 7 2 . However, we will only retain 
first-order terms in k. The latter will be treated on the same footing as 7 2 , so that terms 
in, e.g., 7K will be discarded. 

Considering that quantities of the form — I u )/I p \ are our smallest parameters, we 
will also take those to be of order 7 2 . Accordingly, the quantities T M defined in fl2.5j) will be 
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approximated as T\ ~ A2 — A3 whenever multiplied with 7 or k, and similarly for I^, T3 
with cyclic permutation of the indices. 

In the expansions of the waveforms, linear combinations of A M and T M defined in (|2.5|) 
will appear, multiplied with powers of 7 and k. With the approximations above, these can 
be written as linear combinations of I 3 — (Ii + 1 2 )/2 and 1% — 1\, the first of which measures 
oblateness. Note that we have 

I 2 - h ~ 16« [/ 3 - (h + h)/2], (4.13) 

the correction being of order k 2 . Also, the non-axisymmetry parameter 5 defined in (j3.1b|) 
can be expressed in terms of k as 



1/2 



h-h 

~ 1 - (1 - 16k) 1/2 

~ 8k, (4.14) 



again up to terms in k 2 . 



We are now ready to calculate the spectrum at second order in the wobble angle. 



D. The spectrum at second order 

The expressions for R ifl and f2 M given in Jj| can be expanded to second order in 7 and to 
first order in k. By substituting the results into (|2.8J) and truncating appropriately we find 
the desired second-order expansion of the waveforms. 

In the expansions, the time dependence of the waveforms is given entirely in terms of 
linear combinations of products of trigonometric functions, with frequencies that are integer 
linear combinations of the basic frequencies Q rot and n prec . These trivially allow for a Fourier 
expansion, and after collecting the prefactors of the simple sines and cosines, we find that 
the second-order waveforms h+ and are of the form 

1 

= ^[< fc cos(^)+< fc cos(^t)], 

fc=0 

1 

^ = EK* sin i^ikt) + <t sin (ngf)], (4.15) 

fc=0 
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where we defined 

Qg = tt n + 2ktt pmc . (4.16) 

The contributions with frequencies Q 1 and Q 11 are the first-order ones in small wobble angle 
found by Zimmermann The second-order contributions appear as sidelobes to the first 
order lines. 

We find the following non-zero amplitudes. At first order, we retrieve Zimmermann's 
lines in Jsj]: 

4,Q = -^ 2 (l+cos 2 (i))(/ 2 -/i), 

A% = h 2 sin (2i) [I 3 - (h + I 2 )/2) 7 . (4.17) 

To compare these with the amplitudes from (|3.2|) . note that with the approximations we are 
making, Q pTec <C fi rot , so that b 2 ~ Q 2 ot ~ (f2 ro t + fiprec) 2 (indeed, from ()3.3a|) and ()3.3b|) 
it can be seen that the difference between the three is of second order in 7). The reason 
why we find a simpler frequency prefactor stems from the fact that we have used the closed 
expressions (|2.8jl instead of the more general formulae (j2.1j) as the basis of our computation, 
and that the former do not contain any explicit time derivatives. 
At second order in the wobble angle, we find a single extra line: 

<+i = l b " (1 + cos2 «) [ J 3 " ( J i + J 2)/2] 7 2 - (4.18) 
The results for the 'cross' polarization are quite similar: 

4,o = ~\ h2 cos(i)(I 2 - h), 

A" =-b 2 sin (<) [h - (h + J 2 )/2] 7 (4.19) 



at first order in 7, and 



4,+i = ^ 2 cos (i) [J 8 - (J a + / 2 )/2] 7 2 (4.20) 



at second order. 

Henceforth, the second-order line will be referred to as line III. 

Figure 1 gives a schematic representation of the spectrum at second order in the precession 
angle. 
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FIG. 1: A schematic representation of the spectrum at second order in the wobble angle. The lines 

HQ 

I and II found by Zimmermann et al. |2|, LJ arise at first order in the small wobble angle. Line I is 
due to the combination of rotation and non-axisymmetry. Line II has an amplitude proportional 
to the wobble angle parameter 7. The amplitude of line III is proportional to 7 2 ; it appears as a 
sidelobe to line I. 

V. DETECTABILITY AND PHYSICAL CONTENT OF THE SPECTRUM 
A. Detectability 

In subsection llVBl astrophysical estimates were obtained for the small parameters govern- 
ing the problem: the wobble angle 7, the deviation from axisymmetry k, and the oblateness. 
With the expressions for the amplitudes of the spectral lines found in the previous subsec- 
tion, these can now be used to determine the detectability of the lines. For concreteness, 
we will be interested in distances of r ~ 10 kpc, rotation frequencies / ro t ~ 500 Hz, and a 
rigidity parameter f3 ~ 10 -5 . Factors of G and c in the amplitudes are reinstated. 



15 



The ('cross' polarization) amplitude for line I is of the form 

A 1 = -~/f^{h-h) 

G64 / 3 -(/ 1 + / 2 )/2 

= ~ I ^— h /(*) j 3 ( 5 - x ) 

with /(z) some function of order 1 depending on the polarization; in the second line we have 
used (l4~T3l . Using the estimates of subsection |IVB] and setting I 3 = (2/5)(1.4 M Q )(10 6 cm) 2 , 
we find 

^3.3xHT 28 ( r HJ5£ > ) (VV)^ (^-\\ (5.2) 



500 Hz/ R V10 yr / V.10 3 yr 
For line II we have (under the same assumptions): 

c 4 r i 3 

-26 ( 10 k P<A f Aot Y f P 



4.6x10-1^)1^ ^|7, (5.3) 



where <?(i) is again a function of order 1. Finally, for line III: 



F - 9.4 x 10- 26 ( 'H^£ ) ( ( J- ) 7 2 . (5.4) 

V r y V5oo Hzy Vio- 5 / 

In Jones and Andersson , detectability of spectral lines from precessing neutron stars 
was discussed using estimates of the sensitivity of interferometric detectors that differ from 
the ones we will use. Their paper was written before the more recent analysis of the expected 



capabilities of LIGO II by Fritschel 



14L Here we will assume a LIGO II dotoctor. 



optimized for isolated neutron star observations by adding a signal recycling mirror at the 
output, the position of which can be changed macroscopically. For frequencies between 500 
and 1000 Hz, this leads to a strain sensitivity of 10" 24 Hz" 1/2 Q,Q- A one- year integration 
time would then allow for the detection of signals with amplitudes as low as 

h ■ ~ 10~ 24 ( - ^ ' 

mm " U Ul5 x 10 7 J 

~ 2 x 10" 28 . (5.5) 

From (|5.2|) . we see that line I should be detectable for a spin-down age much less than 10 3 
years. (Also recall that our estimate for the magnitude of k should be considered pessimistic.) 
Eq. (15. 3J) indicates that to observe line II, one would need 7 > 4.3 x 10" 3 . This is to be 
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compared with the physical upper limit implied by eq. (j4.11|) . which is 7 max ~ 1.8 x 1CT 2 
for a breaking strain of Ubroak ~ 10~ 3 . 

Finally, we come to line III. Eq. (|5.4j) shows that one needs 

7 > 4.6 x 1(T 2 (5.6) 

for the second-order line to be detectable. Now, recall that there is much uncertainty about 
the value of the breaking strain of the crust 

Q- e.g., iVeak ~ 10 2 , then for the frequency 
of /rot ~ 500 Hz we are considering, the maximum possible wobble angle would be 0.18 rad. 
Since we have been assuming a distance of r ~ 10 kpc, this means that line III may well be 
observable even for sources as distant as the galactic center. 



B. Extracting physical information 

We now investigate what information about a neutron star's shape and motion can be 
deduced from the spectrum computed in the previous section. 

If the neutron star is not precessing, then at most one line will be present in the spectrum, 
namely the first-order line I associated with non-axisymmetry. We will not consider this case 
any further. 

Even if the non-axisymmetry of the neutron star is insignificant, two lines are present; 
both lines II and III are generated purely by the precessing motion. The inclination angle 
i can always be determined by comparing the two polarizations for two different lines. 
Introducing non-axisymmetry, line I appears, which is separated from line III by two times 
the precession frequency. Let us assume that line III has enough power to be detectable. 
We can then distinguish between two cases. 

(a) The precession period is large compared to the integration time. In that case 
we see the precession line II at a frequency f2 rot + Q prec ~ f2 rot . Even though 
lines I and III may not be resolvable, the first-order non-axisymmetry line I gets 
modulated by the second-order precession line III. When studying data, one 
would know what function to fit the time dependence of this modulation by. 
The unknowns are the precession frequency, and the amplitudes of lines I and 
III. Assuming the fit would allow one to estimate these quantities with some 
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degree of accuracy, the rest of the discussion for this case is identical to that of 
case (b) below. 

(b) If lines I and III can be resolved, three distinct lines are seen. The parameters 
encoding the physical information about the neutron star can then be found as 
follows. The ratio of the amplitudes of lines III and II is proportional to 7, the 
proportionality factor being a function of the inclination angle i. Next, one can 
look at the ratio of the amplitudes of lines I and II. This is proportional to re/7, 
where the proportionality factor is again a function of i; knowing 7, we obtain 
a value for k. The frequencies at which the lines occur allow us to determine 
fl prec and f2 rot . Now, the expression ()3.3a|) for the precession frequency can be 
written as 

"prec — „ n T "rot- \°- 1 J 

2K(m) I 3 

Using m = I6K7 2 , from Q rot and fi prec we can then compute the oblateness 
parameter 

I 3 - (I x + J 2 ) 12 , s 

V . (5.8) 

^3 

It is important to note that, to extract the physical information one is interested in 
(precession angle, deviation from axisymmetry, and oblateness), the first-order spectrum 
would not suffice. 

We end this section with a brief comment concerning the distance to the neutron star 
emitting the radiation. Consider the amplitude of, e.g., line I: 

A' = -I°P*m h - {I \ +h)l \ (5-9) 
c 4 r 1 3 

where f(i) is again some function of the inclination angle of order one. As we have just 
seen, knowledge of the amplitudes and frequencies of the three spectral lines at second 
order suffices to infer values for i, k, [J 3 — (Ji + I 2 )/2]/I 3 , and of course b 2 ~ f2^ ot . From 
the expression above, it is then clear that by inserting an educated guess for the principal 
moment of inertia I3 (as we have been doing to assess the detectability of the lines) , we can 
obtain a rough estimate of the distance to the source. 
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VI. SUMMARY AND CONCLUSIONS 



We derived the quadrupole gravitational wave spectrum of a rigidly rotating and freely 
precessing neutron star up to second order in the precession angle. 

Before doing so, we reviewed Zimmermann's first-order expansion of the waveforms j^] . 
We then introduced new 'small' parameters which are more directly related to the physical 
quantities of interest, with a view on facilitating the second order expansion and its interpre- 
tation, and also clarifying the earlier work. These parameters correspond to the oblateness, 
the deviation from axisymmetry, and the precession angle of the neutron star. 

Next, we made careful astrophysical estimates of physical quantities characterizing neu- 
tron stars. This led to estimates of the 'small' parameters, which justified the use of the 
wobble angle as the primary expansion parameter for a higher-order expansion of the wave- 
forms, putting the other small parameters on the same footing as the square of the wobble 
angle. In addition, the estimates explained why the truncations implicitly made by Zim- 
mermann were physically reasonable. 

We then used the new small parameters to set up an expansion of the waveforms to 
second order in precession angle. The approach taken differed somewhat from that of Zim- 
mermann's in that we used as a starting point the closed expressions for the specific case 
of quadrupole radiation from rigidly rotating and precessing bodies instead of the general 
formulae, truncating directly in terms of parameters with a clear physical interpretation. 
Our method does not imply a significant departure from Zimmermann's approach, but it 
does make the calculations and the end result more transparent. 

At first order, the two spectral lines predicted by Zimmermann were retrieved: line I, 
which is associated with non-axisymmetry, and line II which results from the precessing 
motion. At second order, a single additional line was found (line III). It appears as a 
'sidelobe' to line I and is separated from it by twice the precession frequency. 

On the basis of our estimates for the small parameters, we evaluated the detectability 
of the three spectral lines with a LIGO II detector optimized for isolated neutron star 
observations Q, Q], assuming a one-year integration time. For a rotation frequency of 
500 Hz and a distance of 10 kpc, together with some conservative assumptions concerning 
neutron star structure, we found that: 

• Line I should be detectable if the neutron star's spin-down age is much less than 10 3 
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years; 

• Observability of line II depends on the crustal breaking strain, which is subject to 
much uncertainty. However, the line can be observable even if the breaking strain is 
near the lower limit of the range proposed by Ruderman jsij ] : 

• If the breaking strain is near the upper end of this range, line III can also be detectable. 

Hence, it may well be possible to observe the entire second-order spectrum of neutron stars 
as far away as the galactic center. 

Finally, we investigated what gravitational wave observations might tell us about the 
neutron star emitting the radiation. The two first-order spectral lines found by Zimmermann 
would not allow one to infer much about the intrinsic physical properties of the source beyond 
the rotation and precession frequencies. To separate the precession angle, the deviation from 
axisymmetry and the oblateness from g iVGIl dcltcl, clt lficLSt one additional line is needed. Our 
results provide the missing information: if the second order line can be seen, then it does 
become possible to infer these quantities. 

The direct measurement of these characteristics would permit an evaluation of the models 
or neutron star structure and evolution that led to our estimates of the small parameters 



iy y y 
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2l| , in a way that would be impossible by conventional astronomical 
means. The detailed implications for neutron star modeling will be investigated elsewhere. 

Given our assessment of the observability of the spectral lines we found, it may be of 
interest to extend our results by constructing still higher-order expansions. The scheme 
presented here can easily be used to compute any number of discrete contributions to the 
quadrupole gravitational wave spectrum of freely precessing neutron stars, modeled as rigid 
bodies. 
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